Distribution of the S-matrix in chaotic microwave cavities with direct processes and 

absorption 
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We quantify the presence of direct processes in the S-matrix of chaotic microwave cavities with 
absorption in the one-channel case. To this end the full distribution Ps{S) of the S-matrix, i.e. S = 
y/Re^^ , is studied in cavities with time-reversal symmetry for different antenna coupling strengths 
Ta or direct processes. The experimental results are compared with random-matrix calculations and 
with numerical simulations based on the Heidelberg approach including absorption. The theoretical 
result is a generalization of the Poisson kernel. The experimental and the numerical distributions 
are in excellent agreement with random-matrix predictions for all cases. 
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Random-matrix theory has been successfully applied 
to many different scattering systems in several branches 
of physics ranging from quantum mechanics and meso- 
scopic physics to sound and microwave systems Q, |^ |^ . 
The only precondition is the existence of chaotic ray dy- 
namics of the system. Nevertheless, the connection be- 
tween the statistical properties of scattering and the un- 
derlying chaos is not straightforward since scattering may 
involve two time scales, namely a prompt and a delayed 
response. Prompt or direct processes are those in which 
the waves pass through the interaction region without a 
significant delay. In the equilibrated or delayed processes 
the waves suffer several reflections inside the interaction 
region. The delayed processes are usually studied with 
techniques of random matrix theory, whereas the direct 
processes are described in terms of the average of the 
S'-matrix. 

Startingfrom the pioneering work of Lopez, Mello and 
Seligman |5| , there are several theoretical works address- 
ing the statistical distributions of the S'-matrix with im- 
perfect couphng or direct processes |R R 1^ Holllllll^ . 
This distribution is known in literature as the Poisson 
kernel. In the one-channel case with no absorption it is 
possible to parameterize the S-matrix as S{E) — e*^*^^-*, 
and the Poisson kernel reads 



p{0) 



1 i~i(^)r 
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(1) 



where (S) is the ensemble (or energy) average of S{E). 
For (S) = the distribution of the phase 0{E) is uniform 
between and 27r, i.e., S{E) is uniformly distributed on 
the unitary circle, in agreement with the circular ensem- 
bles of random-matrix theory. Eq. means that the 
S-matrix distribution of a system including direct pro- 
cesses is fixed by the average (S) exclusively. 

Following Brouwer ^3|> the Poisson kernel can be in- 
terpreted as follows: For a chaotic system with an at- 
tached waveguide with ideal coupling, the distribution of 



the S-matrix is uniform. If the coupling in the waveg- 
uide becomes non-ideal then the new S-matrix is dis- 
tributed according to the Poisson kernel. Thus, the Ja- 
cobian of the transformation between S-matrices of the 
system with ideal coupling to the system with non-ideal 
couphng yields the Poisson kernel, Eq. (QJ. 

Unfortunately the microwave experiments cannot be 
directly compared with Poisson's kernel Eq. due to 
losses or absorption. When absorption is present, S is 
a subunitary matrix. For the one-channel case the S- 
matrix can be parameterized as 



S: 



Re 



iB 



(2) 



where R is the reflection coefficient. The coupling be- 
tween the scattering channels and the interior region can 
be quantified by the transmission coefficient Ta of a bar- 
rier describing the direct processes: 



Ta = l-\{S)\' 



(3) 



The subindex a will be used to denote the antenna cou- 
pling. For perfect antenna coupling [Ta — 1) there are 
some exact results. The phase 9 is uniformly distributed 
between and 2n as before. The distribution PRfi{R) of 
R is known in the cases of strong (7 1) and weak 
(7^1) absorption for systems with and without 
time reversal symmetry (/? = 1,2). Throughout this let- 
ter the subindex "0" refers to perfect coupling. In case of 
systems without time reversal symmetry and a single per- 
fectly coupled channel the distribution Pr.o{R) was cal- 
culated for any absorption strength 7 by Beenakker and 
Brouwer as well as for two perfectly coupled chan- 
nels in presence of time-reversal symmetry. For systems 
with imperfect coupling and absorption the distributions 
of the proper time delays and reflection eigenvalues are 
only known for (3 — 2 |1!tI |. More recently, a general dis- 
tribution of the reflection eigenvalues has been obtained 
for a large number of propagating channels, independent 
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FIG. 1: Sketch of the model for scattering with direct pro- 
cesses and absorption. In (a) So describes the scattering of 
a billiard with absorption 7 but perfect coupling to the lead 
{Ta = 1). In (b) we associate the S-matrix Sa given by Eq. 
to the barrier that describes the non-ideal coupling. The re- 
sulting S-matrix of the system with absorption 7 and coupling 
Ta can be written in terms of 5*0 and Sa- 



of time-reversal symmetry 16]. Other quantities, such as 
probabihties of no return, distributions of Wigner time 
delay have also been obtained for systems with absorp- 
tion [HEl. 

From the experimental point of view the effect of the 
absorption has been studied on the ^-matrix correla- 
tion function 0, l^^l, the reflection coefficient 0, Il9| . 
the cross-correlation functions of the ^-matrix j2l| . the 
transmission coefficient |2^ and very recently on the 
impedance matrix |2^ . 

In this paper measurements are presented which pro- 
vide clear evidence of direct processes in chaotic scatter- 
ing with absorption. In particular the distribution Pe{&) 
of the phase of the subunitary S'-matrix given in Eq. ^ 
will be obtained. This corresponds to a generalization 
of Poisson's kernel, Eq. including absorption. Also 
the full experimental distribution Ps{S) of the S'-matrix 
in presence of both absorption and direct processes for 
the one-channel case is given. Both distribution are re- 
lated through (6*) = PsiS)dR. We compare our 
measurements with numerical simulations based on the 
Heidelberg approach and with our random-matrix ana- 
lytical calculations. 

To derive the distribution PsiS) of the S'-matrix in 
the presence of direct processes and absorption for the 
one-channel case we consider the systems of Fig. ^ In 
Fig. n^a) a system with absorption strength 7 and per- 
fect couphng (ra=l) is shown. Let us denote by 5*0 — 
^/Roe^^° the S-matrix describing the scattering in that 
case. We assume that the distribution Ps,o{So) of Sq 
is known. Actually it is sufficient to know the distri- 
bution of the reflection coefficient PRft{Ro) under the 
assumption that Rq and are still uncorrelated when 
absorption is present. As we mention above, the phase is 
uniformly distributed in the case of perfect coupling and 
no absorption, i.e., 



Ps,o{So) 



1 

2^ 



PR.oiRo)- 



(4) 



Just as for systems without absorption, the S-matrix of 
the system with direct processes and absorption given 
by Eq. Q can be written in terms of the S-matrix of 
the system with absorption but without direct processes. 
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FIG. 2: The phase 6 of the S'-matrix for different coupling 
and absorption regimes: (a) 7 — 0.56, Ta = 0.116 (weak 
absorption- weak coupling), (b) 7 = 2.42, Ta — 0.754, (c) 
7 = 8.40, Ta = 0.989, and (d) 7 = 48.00, Ta = 0.998 (strong 
absorption- nearly perfect coupling). The corresponding Ar- 
gand diagrams are shown on the right hand side. 



Sq, and the S'-matrix of the barrier Sa that describes the 
non-ideal couphng of the antenna (see Fig. ^b)). As a 
model, we take the S-matrix of the barrier to be 



Sq 



(5) 



where Ta is the transmission coefficient of the coupling 
barrier. The combination rule of S-matrices gives the 
following relation between S and Sq 



So(S) 



s-{s) 

1-{S)S' 



(6) 



where {S) = VI - Ta. Then the distribution Ps{S) for 
the system including direct processes is Ps,o(So) times 
the Jacobian of the transformation The result is 



Ps{S) 



d{Ro,0o) 



9(i?, e) 

1 - {s? 
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PsASq) 



:^Pr.ARo), (7) 



where Rq — |Sop. The distribution Pb..o{Ro) is known 
for several cases [Ulli: 

2SMii^[A + S(l-i?o)]; /3 = 2, 
PrARo) - { Ce-"/(i-^o)/(i _ i?o)2+/3/2. ^ « 1 

7> 1, 
(8) 
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FIG. 3; Distribution of the S'-matrix for the same ranges of couphng and absorption of Fig. |21 The upper, central and lower 
rows correspond to experiment, numerics and theory, respectively. Note the change of scales for R and 9. 



where a = 7/3/2, A ^ a{e" - 1), B = {1 + a - e"), and 
C = ai+'5/2/r(l + P/2). In the first case of Eq. © the 
absorption can take any value, whereas in the other cases 
f3 can take any value. As our experimental results are for 
a time reversal system we need the distribution for (3=1. 
An approximate distribution for (3=1 was given in Ref. 
|l5l| . that is valid, however, for intermediate and strong 
absorption only. The following interpolation formula 



Pr.o{Ro) — C, 



2-1 



B{1 - R) 



13/2 



(1 - i?)2+/3/2 

(9) 

with Cp = {AT{1 + (3/2,a)/a^ + Be~"/a)-\ where 
r(x, a) = t^^^e^^dt is the upper incomplete Gamma- 
function. It satisfies all cases of Eq. ((SJ. Notice that the 
theory presented here is rigorous for the imperfect cou- 
pling or direct processes, only the ansatz given by Eq. @ 
for the absorption is heuristic. 

The experimental setup was described in Ref. 0| . The 
following two systems were investigated: i) a half Sinai 
billiard and ii) a half Sinai billiard with a microwave ab- 
sorber attached to one side. To improve statistics, in 
both cases the semicircle was moved along the wall in 
steps of 5 mm to obtain more than 50 measurements. The 
complex S'-matrix was measured with a vector network 
analyzer. We investigated four different regimes also con- 
sidered in Ref. ranging from weak absorption and weak 
couphng (7 = 0.56, Ta = 0.116) to strong absorption and 
nearly perfect coupling (7 = 48, Ta = 0.998). The one- 
channel case is realized due to the fact that only a single 
antenna is attached with a radius much smaller than the 



wavelength. 

In Fig. 121 we plot the phase 6* of S* as a function of 
the frequency and the corresponding Argand diagrams 
of the scattering matrix S for different coupling and ab- 
sorption regimes. The experimental value of Ta is ob- 
tained directly from the mean value of the S'-matrix by 
Eq. The measured value of (i?) then fixes the ab- 
sorption strength 7 which decreases monotonically as a 
function of {R) and viceversa. 

To test the theoretical and experimental results we 
made numerical simulations based on random-matrix 
theory using the Heidelberg approach. To obtain the nu- 
merical distributions we used an ensemble of S'-matrices 
constructed according to 

S{E) = I - 2mW\E - H + iT:WW'^)-'^W, (10) 

where H is taken from the Gaussian orthogonal ensemble 
and W gives the coupling between the resonant modes of 
the cavity and the channel states, and E is related to the 
electromagnetic frequency v by Weyl's formula. More 
details of the numerical method can be found in Ref. Q ■ 
In Fig. 13 we plot the experimental distributions of the 
S'-matrix for the same cases as in Fig. |21 The numerical 
distributions are also included in order to test the theo- 
retical ones based on Eqs. Q, 0, and J^. An excellent 
agreement between theory, numerics and experiment is 
found. Finally, in Fig. 01 we show the distribution Pe{d) 
for the different regimes of absorption and antenna cou- 
pling. The numerical simulation is included as well as the 
theoretical results obtained by integrating Eq. numer- 
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FIG. 4: Experimental distribution Pe,T{d) (histogram) for 
the same regimes of coupling and absorption as given in Fig.|21 
Additionally the numerical (crosses) and the theoretical (solid 
line) is plotted as well. In (d) the dotted curves correspond 
to the coupling strengths ra=0.996 and 1, respectively. 



ically. Again an excellent agreement is found for Fig. ^ 
(a)-(c). The deviations in the case of strong absorption 
and nearly perfect coupling, Tq=0.998, (Fig. ^d)) are 
due to the fact that the distribution Pe{0) depends very 
sensitive on the coupling parameter Ta ■ Within the inves- 
tigated frequency range the coupling Ta changes, which 
is not taken into account in the theory. We additionally 
plotted the theoretical distributions for Ta =0.996 and 1. 

In conclusion, we have shown experimental evidence 
of direct processes in the scattering of chaotic microwave 
cavities in presence of absorption. This was done by i) 
obtaining the experimental distribution of the ^-matrix 
for various regimes of absorption and antenna coupling 
in the one-channel case, ii) obtaining the same distribu- 
tions numerically using the Heidelberg approach, and iii) 
presenting a theory that generalizes Poisson's kernel to 
include absorption. Excellent agreement between theory, 
numerics and experiment was obtained in all regimes of 
absorption and antenna coupling. Pr.oI^o) as given by 
Eq. I^is not rigorously derived for (3—1. As soon as 
a rigorously derived result becomes available it can be 



plugged into Eq. j?)) to get the exact distribution of the 
S matrix for any absorption and antenna coupling. 



This work was supported by DGAPA-UNAM, 
Conacyt-Mexico, and the DFG (Germany). MMM 
thanks Centro Internacional de Ciencias A. C. and Cen- 
tro de Ciencias Fi'sicas for kind hospitality. We thank 
C. Lewenkopf, P. Mello and J. Flores for useful com- 
ments. We are especially indebted to D. Savin and Y. 
Fyodorov for suggesting this improved version of Eq. |S1 
using time delay results from Refs. [ill Il7||. 



[1 
[2: 

[s: 
[4: 

[5] 
[6] 

[7 

[8 

[9: 

[10: 
[11 

[12 

[13 

[14 

[15 

[16 

[17 
[18 

[19 

[20 

[21 

[22; 

[23; 



C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 
H.-J. Stockmann, Quantum Chaos, (Cambridge Univer- 
sity Press, Cambridge, 1999). 

P. A. Mello and N. Kumar, Quantum Transport in Meso- 
scopic Systems, (Oxford University Press, Oxford 2004). 
R. A. Mendez-Sanchez, U. Kuhl, M. Barth, C. H. 
Lewenkopf and H.-J. Stockmann, Phys. Rev. Lett. 91, 
174102 (2003). 

G. Lopez, P. A. Mello and T. H. Seligman, Z. Phys. A 
302, 351 (1981). 

T. H. Seligman in Proceedings Cah ELAF 1982 Stochas- 
tic Processes Applied to Physics and other Related Top- 
ics , eds. A. Rueda and S. Moore World Scientific, Sin- 
gapore 1983. 

P. A. Mello, P. Pereyra and T. H. Seligman, Ann. Phys. 
161, 254 (1985). 

W. A. Friedman and P. A. Mello, Ann. Phys. 161, 276 
(1985). 

E. Doron and U. Smilanzky, Nucl. Phys. A 545, 455 
(1992). 

P. W. Brouwer, Phys. Rev. B 51, 16878 (1995). 

D. V. Savin and Y. V. Fyodorov and H.-J. Sommers, 
Phys. Rev E 63, 035202 (2001). 

G. B. Akguc and L. E. Reichl, Phys. Rev. E 67, 46202 
(2003). 

E. Kogan, P. A. Mello, and He Liqun, Phys. Rev. E 61, 
R17 (2000). 

C. W. J. Beenakker and P. W. Brouwer, Physica E 9, 
463 (2001). 

D. V. Savin and H.-J. Sommers, Phys. Rev E 68, 36211 

(2003) . 

D. V. Savin and H.-J. Sommers, Phys. Rev E 69, 35201 

(2004) . 

Y. V. Fyodorov, JETP Lett. 78, 250 (2003). 

Y. V. Fyodorov and A. Ossipov, Phys. Rev. Lett. 92, 

084103 (2004). 

E. Doron, U. Smilansky and A. Frenkel, Phys. Rev. Lett. 
65, 3072 (1990). 

C. H. Lewenkopf, A. Miiller and E. Doron, Phys. Rev. A 
45, 2635 (1992). 

R. Schafer, T. Gorin, T. H. Seligman, and H.-J. 
Stockmann, J. Phys. A 36, 3289 (2003). 

H. Schanze, E. R. P. Alves, C. H. Lewenkopf, and H.-J. 
Stockmann, Phys. Rev. E 64, 065201 (R) (2001). 

S. Hemma dy, X. Zheng, E. Ott, T. M. Antonsen and S. 
M. Anlage, [coiid-mat /0403225 



